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We study the frequency-synchronization of randomly coupled oscillators. By analyzing the con- 
tinuum limit, we obtain the sufficient condition for the mean-field type synchronization. We es- 



pecially find that the critical coupling constant K becomes in the random scale free network, 
P(k) oc k 7 , if 2 < 7 < 3. Numerical simulations in finite networks are consistent with these 



analysis. 
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I. INTRODUCTION 

Recently, it has become clear that the complex network plays an important role in many 
natural and artificial systems, such as neural network, metabolic systems, power supply 
system, Internet, and so on |l|, yj . ^ n particular, we have recognized that many networks have 
scale-free topology; the distribution of the degree obeys the power law, P(k) ~ k" 1 . The 
study of the scale- free network now attracts the interests of many researchers in mathematics, 
physics, engineering and biology. 

The dynamics in the network systems is one of the important themes of the investigation 
of complex network. In this paper, we study the synchronization of the random network 
of oscillators. The phase synchronization in complex network has been studied by several 
authors 3|, \$, while the frequency synchronization has not been studied so much. One of 
the important work on this problem was made by Watts |5j|. He suggested from numeri- 
cal simulation that mean-field type synchronization occurs in small-world network such as 
Watts-Strogatz model. His study was followed by the work of Hong et al, in which phase- 
diagram and critical exponent are numerically studied in detail 6]. These works showed that 
mean-field type synchronization, that Kuramoto observed in globally-coupled oscillators |7[, 
appears also in the small-world network. However, such a study in scale-free network has 
not been performed yet. 

In this paper, we analytically study the frequency synchronization in the random network 
of oscillators. By analyzing the continuum limit of this model, we obtain the sufficient 
condition for the synchronization. Our result shows that in the scale-free random network 
the threshold for synchronization is absent if 2 < 7 < 3. We also carry out the numerical 
simulations and the results are consistent with this analysis. 

This paper is constructed as follows: the next section describes the model of oscillator 
network and derive the continuum limit equation. Section IIHI is devoted to derive the 
sufficient condition for the synchronization from the continuum limit equation. We show 
that the order parameter is different from the one used in previous works, and especially 
we conclude that the threshold for the synchronization disappears in the random scale-free 
network. These results are consistent with the results of the numerical simulations, which 
are described in section HVl In the last section we make a summary of this paper and discuss 
the relation to the other properties of the scale-free network. 



II. OSCILLATOR NETWORK MODEL AND ITS CONTINUUM LIMIT 

First we describe the model we study in this paper. We study the network with N-nodes. 
At each nodes, there exists a oscillator and the phase of the oscillator 9{ is developed as 

~^=^i + K^2 Oij sm(9j - 6i). (1) 

where K is the coupling constant, aij is 1 if the nodes % and j are connected, and otherwise. 
Ui is a random number, whose distribution is given by the function N(u>). 

For the analytic study, it is convenient to use the continuum limit equation. We define 
P(k) as the distribution of nodes with degree k, and p(k,u;t,9) the density of oscillators 
with phase 9 at time t, for given uj and k. We assume that p(k, u; t, 9) is normalized as 

/ p(k,u;t,8)d6=l. (2) 

Jo 

For simplicity, we assume N(u) = N(—u>). Under this assumption, we suppose that the 
collective oscillation corresponds to the stable solution, -£ — 0, in this model. 

Now we construct the continuum limit equation for the network of oscillators. The 
evolution of p is determined by the continuity equation dp/dt = —d(pv)/d9, where v 
is defined by the continuum limit of the r.h.s of eq.(|Tj). Because one randomly se- 
lected edge connects to the node of degree k, frequency u, phase 9 with the probability 
kP(k)N(uj)p(k,uj;t,6)/ J dkkP(k), p(k,u;t,9) obeys the equation 



dp(k,u;t,9) d 

di ~~d9 

In the next section, we study the mean-field solution of this equation. 



„J KkJdcu'Jdk'Jd9'N(uj')P(k')k'p(k',u;';t,9')sm(9-9') 

p(fc ' " ; *' e) r + Tdypww> 

(3) 



III. MEAN-FIELD ANALYSIS OF RANDOM OSCILLATOR NETWORK 

In this section, we study the sufficient condition for the synchronization using eq.©- 
First we introduce order parameter (r, ip) as 

re** = fdufdk f d9N{u)P{k)kp{k, u; t, 9)e ld / f dkP{k)k. (4) 

This order parameter is different from the one which is used in previous work in small-world 
modelJ5|, |6(. In the previous works, J2i et9% /N is used for the mean-field, while our order 



parameter corresponds to J2i k%e t0i /YA h, where ki is the degree of the node i. However, 
from eq.Q it seems natural to use eq.(j3J) as the mean-field value the in random network. 
Here we note that < r < 1. 

Inserting eq.fjlj) into eq.Q, we get 

dp ^^ 9) = ~{p{k, w; t, 9){u + Kkr sin(^ - 0))). (5) 

The time-independent solution of p is then 

^(p(k, w; t, 9)(u + Kkr sin(</> - 9))) = 0. (6) 

Without a loss of generality, we can assume ip — 0. Since we want to seek the solution 
which corresponds to the Kuramoto's solution in globally coupled oscillators, we assume the 
solution of this equation as 

f 8(9 - arcsinf-^-)) if #- < 1 
p{k,u-9) = \ \ {Kkr)) « k '- (7) 

1 — eV ■ ai otherwise, 

^ \u— Kkr smti\ 

where C(k,u) is the normalization factor. Here we note that p depends on both K and 
k. This equation means that Kk corresponds to the coupling between mean-field and the 
oscillator. Inserting eq.(jZJ) into eq.(jlj), we get the equation for r, 

r = fdufdkf d9N{uj)kP{k)p{k, u; 9)e ie / f dkkP{k). (8) 

To calculate this integral, first we divide the integral over u. 

//■ /■ f f ( fKkr r—Kkr poo \ 

du dk d9N(u)kP(k)p(k, u; 9)e ie = dk d9 / du + / duo + du\ 

x N(uj)kP(k)p(k,uj;9)e ie (9) 

The contribution from the integral at u> < —Kkr and uj > Kkr is if N(u>) = N(—oo), 
because 

—Kkr poo poo I I 

du+ du)N(u)p(k,u;9)e i6 = / N{u>)e»C(k,u)( . + . ) 

-oo ./.ftrfcr ./iffer c<j — A At sin u; + Arer sin 

(10) 
The integral of the r.h.s of eq.(fTUj) over 9 equals to 0. Therefore eq.® is equivalent to 

/pKkr fjj p 

dk N(u)kP(k) exp(i arcsinf — — ■))/ dkkP(k). (11) 

J-Kkr Kkr J 



If we assume arcsin(^-) is between [— vr/2, vr/2], we get 



Kkr / ,,, \2 



dk / duN(u)kP(k)Jl 



( " X 

\Kkr, 



I dk f duo'kP(k)N(Kkruo')^l - uo' 2 x Kkr 
Krfdkk 2 P(k)f du'N(Kkru')Vl - to' 2 . (12) 



If r ^ 0, we get 



fdkkP{k) = K fdkk 2 P{k) f du'N{Kkru')Vl-oo' 2 (13) 

The l.h.s of this equation is independent of r and we define the r.h.s of this equation as f(r). 
At r = 1, /(r) is not larger than / dkkP(k), because 

f dkk 2 P{k) f du'N(Kkr)Vl-oo' 2 < f dkk 2 P{k) f duj'N(Kkru') 

/I r°° 
dkk 2 P{k)—— / duo"N{uo") 

JdkkPjk) 

Kr ' { } 

here we use the relation f^duN^u) = 1. Therefore the sufficient condition that eq. (|12|) 
have solution at < r < 1 is that f(r) > J dkkP(k) at r = 0, 

KN(0)7rJdkk 2 P(k) 

2JdkkP(k) [ ' 

This is the sufficient condition for the synchronization in random network of oscillators. The 
most impressive point of this equation is that in the random scale-free network, P(k) oc A;" 7 , 
this condition is satisfied for any fC>0if2<7<3, because / dkk 2 P{k)/ J dkkP(k) 
diverges. Therefore we have no threshold for the synchronization in the random scale- 
free network. This seems similar to the absence of the threshold in susceptible-infected- 
susceptible(SIS) model J8j. We will discuss on this similarity in later section. 

In this section, we derive the sufficient condition for the synchronization in random net- 
work of oscillators, using the continuum limit equation. In the next section, we show that 
the analysis above is in good agreement with the results of the numerical simulations. 

IV. NUMERICAL SIMULATION OF SYNCHRONIZATION 

In this section, we show the result of the numerical simulations of the random network 
of oscillators. In the all simulations, we take N(u) as N(u) = 0.5 if —1.0 < u < 1.0, and 




(a) K = 0.15 



(b) K = 0.30 



FIG. 1: (uj,d6/dt) distribution of oscillators in random network, for K=0.15 and K= 0.30. 



otherwise. 

First we show the result on the 1000-node Erdos-Renyi random network model. We choose 
the probability of coupling p = 0.005, which gives / dkkP(k) = 5.0 and / dkk 2 P(k) = 29.7 
on average. In this case, estimated critical K is K c = 0.214. Each simulation is carried out 
100 times. 

In FigQwe plot the relation between uji and ^ after a long time (t = 200) when K = 0.15 
and 0.30. In the case of K = 0.15, d9/dt seems to depend on uj linearly. On the other hand 
at K = 0.30 many oscillators seem to be synchronized at dd/dt = 0. This figure strongly 
suggests that the synchronization occurs between K = 0.15 and 0.30. 

We plot the relation between uj and 9 for K = 0.15 and K = 0.30 in Fig. El We find a 
clear difference between these two cases. In the case of K — 0.30, the distribution of (a>i, 9j) 
is apparently non-uniform, while at K — 0.15 we cannot find any structure. In the case of 
K = 0.30, 9 seems to depend linearly on uj. However, from the previous analysis we suggest 
that 9 depends on both uj and k. To clarify the degree dependence, we plot (a;, 9) for the 
nodes with the degree k = 3, 5 and 7 at K = 0.30 in Fig. El- We also plot aicsm(uj/Kkr) in 
these figures. The average of r is 0.623 in our simulation. From these figures, we find that 
distribution of (uj, 9) seems to be concentrated around a single line. The concentration line 
qualitatively coincides with 9 = aicsm(uj/Kkr). This result suggests that our mean-field 
defined by eq.(JU) is the correct one. 

To estimate the critical coupling K c , we plot K dependence of the average of the order 
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(a) K = 0.15 



(b) K = 0.30 



FIG. 2: (w,0) distribution of oscillators in random network, K=0.15 and K= 0.30. 




(a) k = 3 



(b) fc = 5 



(c) fc = 7 



FIG. 3: (w, 9) distribution of oscillators with degree 3, 5 and 7 in random network, K= 0.30. 



parameter r av . in Fig0] r av . is less than 0.1 and shows weak dependence on K at K < 0.2. 
This non-zero value of r a v. is due to the finite-size effect. On the other hand, at K > 0.2, 
r av . increases rapidly as the interaction increases. This figure suggest that K c is about 0.2, 
which is in agreement with our analysis. Therefore we conclude that all numerical results 
are consistent with our analysis. 

From these simulation, we find that our mean-field theory is applicable to Erdos-Renyi 
model. However, the most impressive suggestion of our analysis is the absence of the thresh- 
old in the random scale-free network. In the following, we show the result of simulation in 
random scale-free network with 7 = 2.5. In Fig. we show the relation between order 
parameter r and coupling constant K for N=500, 1000, 2000 and 4000. At N=500 , r av . 
rapidly increase above K ~ 0.16, which is qualitatively consistent with the K c ~ 0.175 
estimated from eq. (fT2"j) . As the network size increases, r av . at small coupling decreases, 




FIG. 4: Interaction dependence of mean-field parameter rav.- 

which suggests that the finite r av . at small coupling is the finite-size effect. The order pa- 
rameter begin to increases rapidly above K c . We note that by increasing the system size, 
the increase of order parameter begins at smaller coupling. This means that the critical 
coupling K c decreases as the system size increases. We also show K c estimated from eq. (|12|) 
in this figure. The estimated K c qualitatively coincide with the coupling constant at which 
the order parameter increases rapidly. We conclude that our analysis and the results of the 
numerical simulation show a good agreement also in the random scale-free network. These 
results suggest that in the infinite size scale-free network the critical coupling constant K c 
becomes zero, just same as the continuum limit equation. 

To compare the results of the numerical simulation and the analysis more precisely, we 
need a more accurate estimation of K c from the numerical simulation. In the case of the 
globally coupled networks and Watts-Strogatz model, K c is numerically obtained as the point 
at which iV a25 rav. becomes independent of the size of the network [a]. In their analysis, there 
exists an assumption that K c does not depend on the size of the network. On the other hand, 
our analysis and simulation show that K c depends clearly on the size of the network through 
the average of the square of the degrees. Therefore it is impossible to obtain accurate K c 
from finite-size analysis. The exact estimation of K c is a difficult task. 

However, we find that K c derived from eq. (jl2j) seems to have a strong relation to the 
phase transition. We rescale K by K c which is obtained from eq. (j!2j) . and plot the relation 
between iV a25 r av . and K/K c in Fig© In the case of N = 1000, 2000 and 4000, the well- 
defined crossing point exists at K/K c = 1. In the case of N = 500, iV 0,25 rav. a t K/K c = 1 
is a little larger than in the other cases. However, this difference is small, and it seems 
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FIG. 5: Interaction dependence of the mean-field parameter in the random scale-free network for 
N=500, 1000,2000 and 4000. The arrow shows K c estimated from ea. (|12(l . Simulations for each 
parameters are carried out at least 50 realization of the networks. 




FIG. 6: The relation between _/V a25 r a v. and K/K c , where K c is the value estimated from ea. (|12j) . 



that K/K c = 1.0 is the crossing point at large N. This result is similar to the results of 
the finite-size scaling in the globally coupled networks and Watts-Strogatz model. In these 
models, there exists a crossing point at K = K c . On the other hand, our analysis is not 
a precise determination of the critical coupling strength. To avoid the size dependence of 
K c , we rescale K and we have no guarantee that such a rescaling is valid for scale-free 
network model. However, our result strongly suggests that K c obtained from the numerical 
simulation coincides with the result of the analytic solution. 

To conclude this section, we carried out the simulations on Erdos-Renyi model and the 
random scale-free network. All the results of these simulations show a qualitative agreement 
with the analysis in the previous section. 



V. SUMMARY AND DISCUSSION 

In this paper, we study the frequency synchronization of the random oscillator network. 
By analyzing the continuum limit equation, we find that mean-field type synchronization 
occurs in random network model. We obtain the sufficient condition for the synchronization. 
Especially we find that the threshold for the synchronization is absent in scale-free random 
network if 2 < 7 < 3. The results of numerical simulations are in good agreement with this 
analysis. 

One of the most astonishing result in the dynamics of scale-free network is the absence 
of epidemic threshold in SIS model. Our result seems to be similar to the result in SIS 
model, however, there exists a large difference between them. In SIS model the absence of 
epidemic threshold is the result of the divergence of fcnn, the mean degree of the nearest 
neighbor nodes [9(|. On the other hand, in our model the absence of threshold originates 
from the degree- dependence of the coupling between order parameter and oscillators. The 
coupling between the oscillators and the mean-field is proportional to the degree of the 
nodes, as shown in eq.(|5|). and the contribution to the order parameter from the oscillator 
is also proportional to the degree of the node, shown in eq.(|3J). These degree-dependence 
result in the fc 2 -dependence of eq. (|15|) , which leads to the absence of the threshold in random 
scale-free network. Therefore there is a large difference between the absence of threshold in 
SIS model and the synchronization, though these are apparently similar results. To clarify 
this difference, we will need to study the synchronization in the other network models. 
As Eguiluz and Klemm has shown, scale-free network with large clustering coefficient has 
epidemic threshold in SIS modelllOJ, due to the smallness of /cnn- The different behavior of 
the threshold may appear in our oscillator network, because the absence of threshold is not 
caused by the divergence of k nn , but by degree dependence of mean field-oscillator coupling. 
The study of the synchronization in other scale-free network models is a future problem. 
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